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Motivated by recent experiments, we study the interaction corrections to the damping of magne- 
tooscillations in a two-dimensional electron gas (2DEG). We identify leading contributions to the 
interaction-induced damping which are induced by corrections to the effective mass and quantum 
scattering time. The damping factor is calculated for Coulomb and short-range interaction in the 
whole range of temperatures, from the ballistic to the diffusive regime. It is shown that the domi- 
nant effect is that of the renormalization of the effective electron mass due to the interplay of the 
interaction and impurity scattering. The results are relevant to the analysis of experiments on mag- 
netooscillations (in particular, for extracting the value of the effective mass) and are expected to 
be useful for understanding the physics of a high-mobility 2DEG near the apparent metal-insulator 
transition. 

I. INTRODUCTION 

The influence of electron-electron interaction on transport properties of low-dimensional disordered conductors at 
sufficiently low temperatures T remains one of central topics of the condensed matter physics. In a seminal set of 
works (see the reviev^l), Altshuler and Aronov studied the effects of the interplay of interaction and disorder on 
conductivity and tunneling density of states in the diffusive regime characterized by the condition 2ttTt -C 1, where 
r is the transport mean free time (we set fcg = S = 1). Their results were generalized within the framework of 
the renormalization group (RG) by Finkelstein 2 . The last decade has witnessed a renewed increase of activity in 
this field, largely motivated by experiments on an apparent metal-insulator transition in 2D systems. This interest 
was triggered by experiments 3 which showed a "metallic" behavior (decrease of resistivity with lowering T) in high- 
mobility Si structures. Later, qualitatively similar behavior was observed in a variety of high-mobility 2D systems, 
see Refs. 0|H3 for reviews. 

The metallic behavior has been attributed to the effects of the electron-electron interaction in the ballistic tem- 
perature range, ItxTt 3> 1. These effects were originally considered in the framework of the temperature-dependent 
screening 8 . More recently, a systematic theory was developed, taking into account also exchange contributions and 
the effects of both parallel and transverse magnetic fields, and valid in the whole range of T from the diffusive to 
the ballistic regim oV . Another mechanism that can explain the metallic behavior of resistivity in an intermediate 
temperature range in the diffusive regime was studied within the RG framework in Ref . IllL It is applicable to systems 
with more than one valley, such as silicon MOSFETs. 

Despite these successes of the theory, numerous experimental observations remain puzzling and wait for an expla- 
nation. In particular, it was found that the spin susceptibility, proportional to the product mg of the effective mass 
m and the g-factor, is strongly growing when the density approaches the value n c corresponding to the apparent 
transition. This conclusion was drawn on the basis of several experimental methods, including the analysis of beat- 
ing pattern of Shubnikov-de Haas oscillations- 2 '-- 3j - 14 '— , the study of magnetoresistancsi&il in the parallel field, and 
measurement of thermodynamic magnetization^; see recent reviewsi&SS. 

The enhancement of susceptibility with lowering density, interpreted in a number of papers as its divergence at 
ri = n c , has attracted a great deal of attention, since it might be an indicator of some phase transition that the 
system undergoes with a decrease of density. The interpretation of the data has remained, however, controversial. In 
particular, it remained unclear whether the strong increase of spin susceptibility should be attributed to that of m or 
of g. This information is of crucial importance for understanding the nature of the possible transition. 

Several experimental approaches have been used to separate the behavior of the effective mass from that of the 
g-factor. In Refs. 21 22.23 a fit of the resistivity data to the theoretical formulas of Ref. was used to find the 
interaction constant Fq, and thus the g- factor. The accuracy of this procedure is questionable, since the theory of 
Ref. neglects higher Fermi- liquid interaction constants and assumes isotropic impurity scattering. Another approach 
is based on thermodynamic measurements in strong magnetic field24. However, the authors of this work were able 
to measure the effective mass in a very narrow interval of electron concentration only, so that the results are not too 
informative. Also, a strong magnetic field is expected to influence strongly the characteristics of the electron liquid, 
so that the applicability of such measurements to the low-field properties is questionable. So, while most of the above 
measurements seem to indicate that it is the effective mass that is responsible for the strongly enhanced susceptibility, 



2 



an independent verification is clearly needed. 

A well known method for determination of the effective mass is based on the investigation of the temperature 
dependence of Shubnikov-de Haas oscillations (SdHO). It was applied to the present problem in E.efs. H3l25l . However, 
the analysis of the SdHO data is complicated by the fact that both the effective mass and the elastic quantum scattering 
time T q are T-dependent, in view of the combined effect of interaction and disorder. An unambiguous interpretation 
of experimental data requires^ a theoretical information on T-dependence of m and r q . A development of the 
corresponding theory is the aim of the present paper. 

In fact, a recent paper— has made an important step in this direction. Specifically, it was shown in Ref. [2(| that 
the Lifshitz-Kosevich formalism 27 , originally developed for the analysis of magnetooscillations in a 3D Fermi liquid, 
is also applicable in 2D in the regime where the oscillations are exponentially suppressed by temperature smearing or 
disorder. (In the regime of strong oscillations, the Lifshitz-Kosevich formula in 2D should be modified, as was earlier 
shown in Ref.[28|.) Another result of Ref. is that the inelastic electron-electron relaxation does not contribute to the 
damping of magnetooscillations (similarly to the earlier result of Ref. |2^ for the case of electron-phonon scattering) . 

The authors of Ref. then calculated the contribution to the damping induced by the interplay of interaction 
and disorder. Their theoretical treatment of the problem is, however, far from complete. First, they consider only 
diagrams for the self-energy with one impurity-ladder vertex correction to the interaction line and discard diagrams 
with no and with two vertex corrections. Second, they claim that the T-dependence of the oscillation damping rate 
can be equivalently attributed either to the correction to the effective mass, or to the quantum scattering rate (Dingle 
temperature). Furthermore, in the latter case their result for the T-dependence of r q is in contradiction with the 
picture of Friedel oscillations inducing a correction to the relaxation rate, which is linear in T and is governed by 
backscattering^. 

In addition to the above experimental motivation, the development of the theory of interaction effects on magne- 
tooscillations in a disordered 2DEG represents a fundamental theoretical problem. Such a theory should complement 
the recently developed theory of interaction effects in transport of 2D electrons in zero and non-quantizing magnetic 
field o?' 1 ^ . Let us emphasize a peculiar aspect of the present problem. The damping of oscillations is governed by the 
self-energy, which is a single-particle quantity. [Indeed, the relevant diagrams, see Sec III Bl below, are reminiscent 
of those for the tunneling density of states (DOS)]. Generally, the self-energy is not a gauge-invariant object. On 
the other hand, the magnetization and the conductivity (magnetooscillations of which we would like to study) are 
observable (and thus gauge-invariant) quantities. It is well known that the gauge-invariance is of crucial importance 
for interaction-induced corrections; a difference between the results for conductivity and for tunneling DOS in the 
case of Coulomb interaction^ serves as a nice illustration. It is thus a theoretical challenge to see how the gauge 
invariance manifests itself in the magnetooscillation problem. 

The outline of the article is as follows. The section [U] is devoted to presentation of the general formalism. In 
Sec. IIIII we apply it to calculate the interaction-induced contribution to the damping of magnetooscillations in the 
case of short-range interaction. In Sec. lIVI we show how to extract from the above result the corrections to the effective 
mass and the quantum scattering time. We also perform a calculation of the correction to the scattering time based 
on the picture of Friedel oscillations and demonstrate a complete agreement between the two approaches. In Sec.lVlwe 
generalize our results to the case of Coulomb interaction. We also perform there a comparison with the calculation of 
Ref. l26t Section IVII summarizes our findings. Some technical details of our calculations are presented in Appendices. 



II. MAGNETOOSCILLATIONS: GENERAL FORMALISM 



A. Derivation of the formula for a decay of the oscillations 



We begin by calculating the oscillatory part f2 osc of the thermodynamic potential Q. From this quantity one can 
derive the oscillating contribution to the thermodynamic density of states 

r)n B 2 Q 

where /i is the chemical potential, and de Haas-van Alphen oscillations of magnetic susceptibility 

_ u ll osc 

Xosc — , \.&) 

where B is a magnetic field. The main subject of our interest is the exponential damping factor of these magnetooscil- 
lations. For non-interacting electrons, the same exponential damping factor governs the magnitude of the Shubnikov - 
de-Hass oscillations of the conductivity for the case of weak disorder potential in sufficiently weak magnetic field 35,36 , 
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where the self-consistent Born approximation 35 (SCBA) is valid. As we are going to show, the interaction-induced 
correction to the damping factor of the thermodynamic density of states arises due to the renormalization of the ef- 
fective mass and the quantum scattering time. Therefore these T-dependent corrections to the damping factor govern 
the magnitude of the Shubnikov - de-Hass oscillations as well, similarly to the non-interacting case. 

Our starting point is the expression for the thermodynamic potential derived in the paper by Luttinger and Ward^S 



fi = -TTrlnf-G- 1 ) -TTr(GE) + Q'. (3) 

Here 

G(ie n ,muj c ) = [Gq 1 (ie n , muj c ) - T,(ie n , raw c )p' , (4) 

is the dressed Green's function in the Matsubara formalism, ie n = (2n+ l)inT is the Matsubara fermionic energy, lu c 
is the cyclotron frequency, and m is the Landau level index. Further, 

Go(ie n ,mLu c ) 



ie n + Li — (m + l/2)w c 



is the Green's function in the absence of disorder and interaction, and S(ie„,mcj c ) is a self-energy part of Green's 
function which includes all the disorder and interaction effects. 

The trace in Eq. @ implies summation over Matsubara frequencies e n and over Landau levels m. The logarithmic 
term contains all the closed loop diagrams with insertion of self-energy (Fig. 0. The terms — TTr(GE) and Vi' are 
introduced to avoid double-counting of diagrams^^i. The term SI' denotes the sum of all so-called skeleton diagrams 
with all bare Green's functions replaced by dr essed Green's functions (for the recent discussion of Luttinger- Ward 
formalism in 2D Fermi systems see Refs. 1321331) . 

As shown in Ref. the exponential decay of magnetooscillations is described by the Tr ln-term. The oscillatory 
parts of the additional terms, which are introduced to fight overcounting, cancel each other. In order to obtain the 
correction to the thermodynamic potential we need to calculate the self-energy part of the Green's function. 

We decompose the self-energy into two parts: 

E(ie„, mu c ) = Tf dis (ie n , mu> c ) + S cc (i£„, mu c ). (5) 

where £dis(*£n, "^c) denotes the self-energy part due to the scattering on disorder potential with electron-electron 
interaction switched off and S ee (ie n , mw c ) contains all the interaction effects. 

In this paper we assume that disorder potential is <5-correlated, inducing a large-angle scattering of electrons. The 
disorder-induced (noninteracting) part of the self-energy for white-noise disorder and weak magnetic field lo c t <C 1 is 
given by 

hdis(i£n,muj c ) — — . (6) 

ZT 

For stronger magnetic field (i.e. for separated Landau levels), one should employ the SCBA. In this case, both the 
real and imaginary parts of the self-energy depend on ie n in a non-trivial way. In this paper, however, we will address 
only weak magnetic fields when Landau levels overlap. 

>'. 

-TTr ln(— G -1 ) = -T Tr ln(— G^" 1 ) + 

FIG. 1: The logarithmic term £7, Eq. Q, in the thermodynamic potential, Eq. is a sum of closed loop diagrams with 
self-energy insertions E. This term is responsible for magnetooscillations. 

We thus consider the relevant term = — TTrln(— G _1 ) in the thermodynamic potential, 

oc oo 

O = —2v T lu c In [uj c m — fi — ie n + S(ie„, ui c m)] . (7) 

n— — oo m— 
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For overlapping Landau levels, 



UJ C T < 1, 



(8) 



the fc-th harmonics of the magnetooscillations 



_ \ - 2kir 2 n e 
"osc = 2 , A k cos ■ 



k 



eB 



is damped by disorder even at zero temperature via the standard Dingle factor exp(— tt/uj c t). Therefore we will 
consider only the first harmonics of the oscillations, Ax, neglecting all Ak with k > 1 (whose damping is much 
stronger) . 

The oscillatory part of Q is calculated in Appendix A: 



f2 osc ~ 2v y— J A x cos 



2TT 2 n e 
eB 



(9) 



(n e is the electron concentration) with the amplitude of the first harmonics of the oscillatory part of the thermodynamic 
potential given by 



A x = 



4tt 2 T 



£„>0 



exp 



2vr 

UJ* 



1 



2r(l + a ) 



(10) 



Here 5Y<(ie n , £o) is the self-energy part related to the interplay of disorder and interaction. It is analytically continued 
from the points moj c to the whole complex plane £ and taken at £ = £0 where £o [defined in Eqs. (|A.9|) and l|A.10|) 
of Appendix A] is the pole of the Green's function in the presence of disorder. 

The coefficients /?o and «o determine the Fermi-liquid (FL) renormalization of the effective mass in a pure system 
at zero T, 



1 



m = m- 



a 



i + /V 

The effective mass m* in turn governs the expression for the FL-renormalized effective cyclotron frequency, 



eB 



i + A) 



c * L 1 i 

m* 1 + ao 

The coefficient ao is related to the FL renormalization of the Z-factor, 

1 

Z = 



1 + a 



(11) 



(12) 



(13) 



which is given by the residue of the Green's function. 

Depending on the relation between temperature T and the elastic scattering rate 1/r, there are two regimes: 
ballistic, Tr 3> 1, and diffusive, Tr <C 1 (more accurately, the relevant dimensionless parameter is 2itTt). In the 
ballistic regime, it follows from Eq. (JSJ that T ^> 1/r ^> lo c . The diffusive regime can be further split into two 
sub-regimes: normal diffusive (oj c -C T <C 1/r) and ultra-diffusive (T <C lo c -C 1/t). When T ^> cj c , as in the 
ballistic and normal diffusive regimes, only the first Matsubara frequency £o = 7rT in the sum determining A\ is 
relevant, since the contribution of higher Matsubara frequencies are exponentially suppressed. On the other hand, in 
the ultra-diffusive regime T«w c and higher Matsubara frequencies contribute as well. 

In what follows we concentrate on the case T ^> u> c . Under this condition, we get 



Ax 



4n 2 T 



■ exp 



2tt 



{ttT + iZ6Z(iirT,£ )} 





7T 


exp 






UJ*T* 



where we introduced the FL-renormalized scattering time 

T* = T(l 



a ). 



(14) 



(15) 



We note that this renormalization of r is incomplete since it does not include the FL vertex corrections to the impurity 
scattering line. The corresponding contributions is contained in <5£(i7rT, £ ) [T-independent terms in Eqs. (|41|l and 



(a) 




'01 




E 




'00 



'01 
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FIG. 2: Self energy diagrams in the first order in the effective interaction (wavy line). Black triangles denote impurity ladders 
r dressing interaction vertices, Fig. |3] (a) "Simple" self-energies (b) Hikami-box self-energies E^ generated by covering 
each of Ey by an impurity line (dashed) . 
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FIG. 3: Interaction vertex renormalized by the impurity ladder, T(iuk, q)- Dashed line represents scattering on impurity. 



1)8 1( 1 ] and will be addressed in Sections III and IV below. We also note that the inelastic contribution to the self- 
energy oc [(7rT) 2 — e 2 ] sgne„ (see Appendix B) vanishes for e » = ttT , and thus does not affect the damping of the 
magnetooscillations B(T) for T ^§> uj c , in agreement with Refs. I26l29l 

Using the renormalized quantities in Eq. 1(141) we represent A\ in the form 



Ay =4 0) Cn exp[B(T)], 



where 



2n 2 T 



is the standard FL Lifshitz-Kosevich result and 

B(T) 



2TTiZ6Y<(iirT,£ ) 



(16) 



(17) 



(18) 



We thus see that in order to obtain the additional interaction-induced damping factor of magnetooscillations it is 
necessary to evaluate Z 8Yi{mT, £ ). 



B. Self-energy 



We begin by considering the interaction- induced self-energy part Ti ee (ie n , muj c ) in the lowest order in interaction. 
This is sufficient in the case of a weak short-range interaction analyzed in Sec. IIHI below. For the more realistic 
case of the Coulomb interaction (Sec. 0, the relevant higher-order terms can be treated using the random-phase 
approximation (RPA). Higher-order contributions to the T-dependent part of the self-energy, <5£(i7rT, £o)j ar e small 
in the parameter 1/Epr or T/Ep. 

Let us list important elements which are necessary for calculation of the interaction-induced part of the self-energy 
(Fig. |2J). Each contribution to the self-energy has exchange and Hartree parts. We first address the exchange 
contribution (the Hartree terms can be written in a similar way). It contains the angle-averaged Green's function 
covered by the effective interaction line. The corresponding vertices may be dressed by impurity ladders (Fig. |5). 
Notice that the renormalized vertex includes at least one impurity line. Finally, when the interaction line changes the 
signs of Matsubara frequencies at vertices, an additional diagram (we term it a "Hikami-box diagram") with a single 
impurity line covering the whole block is to be included (Fig. [2b) • 
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We split T, cc (ie n , mu c ) into three contributions, corresponding to different possibilities of dressing the two interaction 
vertices by impurities, 

S co = Sqo + 2Soi + Sn. (19) 

Here the subscripts i,j — 0, 1 indicate whether the corresponding vertex is dressed by an impurity ladder. The factor 
2 in front of Sol term reflects two possibilities of dressing one of the interaction vertices. Each of the terms E^ is a 
sum of two contributions, E^- = E"- + E^ , where is the "simple" self-energy (Fig. Et) and E^ is its Hikami-box 
counterpart (Fig. |2p). We note that in Ref. l26l only one out of six diagrams (namely, Egi) was taken into account. 
The expression for Eq in a finite magnetic field reads 

/d 2 q 
- — -2 Jf L ,(qR c )V(iu>k,q) G (ie n - iu k ,mu> c + Luj c ) , (20) 
(27T) 



u> k L 



where 



G{ie n ,muj c ) = - — — — — — (21) 

ie n + jU — (m + l/2)u! c — E(ie„, mu) c ) 

is the Green's function in Landau levels representation, V(iiv k ,q) is the effective interaction, and J? L ,(qR c ) describes 
the bare vertex function connecting Landau levels to and to + L in the quasiclassical limit m,m + I> 1 (for details 
see, e.g., Ref. ^ and references therein). In the expressions for Ej^ and E^ this bare vertex function is multiplied 
by r(iw/ £ ,q) and T 2 (iwfc,q), respectively, where T is the impurity ladder (Fig^J). 

In the limit B — > 0, the corresponding self-energies depend on the momentum p instead of the Landau level index 
to. The interaction vertices are dressed only when the Matsubara energies at the vertices have opposite signs, which 
restricts the summation over u k in £oi an d Sn to the domain E n (u> k — £ n ) > 0: 

Eoo(fe*,&) = -TV I -^F(iw fc ,q) G(fe n -iw fc ,p-q), (22) 
s oi(* £ ™>Cp) = ~T / "T2 v ( iuj k,<l) r(icj fcl q) G(ie n -iu k ,p-q), (23) 



e„(u)*-6 n )>Q' 



E?i(»e„,e P ) = -T £ [A* V (iu k ,tti r 2 ^ fc ,q) G(ie„-iw fc) p-q), (24) 

. W n •/ (27T) 



e n (u)k-s n )>Q' 



where G(ze n ,p) = [ie n + fi — £ p + isgne n /2r — T, ee (ie n , 1 with £ p = p 2 /2m and the vertex correction (FigJSJ 
reads 

r^ fc ,q) = ^ 1 _ . (25) 

^(\ij k \T + l) 2 + {qv F T) 2 -l 

To calculate the damping factor of the oscillations, we need the self-energy taken at the value of £ which is determined 
by the pole of the Green's function, £, = £,o, see Appendix B. According to Eq. (|18|) . the self-energy is further multiplied 
by Z in the damping exponent. This is equivalent to calculating the following integral: 

J 2m 1 + a J J 

~ sgne„^ J dkG{ie n ,k)^{ie n ,i k ). (26) 

We note that the Z-factor drops out in the product ZE. Indeed, the Green's function under the interaction line in 
the self-energy contains the Z-factor in the numerator so that in the numerator of the product ZH we get the factor 
Z 2 . However, the Z-factor is not a gauge invariant quantity and therefore should not appear in the expressions for 
observables, in contrast to the FL-renormalized effective mass. At this point we should take into account the FL 
renormalization of the two interaction vertices in E. Since we are interested in the contribution of relatively slow 
transferred momenta and frequencies giving rise to the T dependence of B(T), q <C kp and tv k <C Ep , we can set 
them to zero when considering the FL vertex renormalization. Then one can apply the Ward identity for the FL- 
interaction dressing of the vertices^, which amounts to multiplying each vertex by a factor 1/Z. These vertex factors 
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in the denominator of B(T) cancel Z 2 in the numerator of B(T). This implies that one can simply discard such 
renormalizations, setting Z = 1 everywhere, when the observable quantities are calculated. The FL renormalization 
then amounts to replacement of the bare band mass m, Fermi velocity v F , and elastic scattering time, r, by the 
renormalized parameters, to*, v* F and r*, respectively. In what follows, we will omit the asterisks, using the notation 
to, vf and t for the renormalized quantities. 

The relevant contributions to the self-energy are calculated in Appendix B. Combining all the terms together, we 
have 



K(w m ,q) = 



d 2 q 

.... - ( 2 -) 2 
[l+r(iw ro ,q)] 2 



V(iu m ,q)K(ico m ,q), 



S(iu m ,q) 



1 - 



W 



where 



and 



S(iu> m ,q) 



1/rf 



? 9 

v F q 



rS 2 (iuj mi q) 



W 2 + v 2 F q 2 , 



1 



S (iuj m ,q) 



W 



1/t 



S Q (iuj m ,q) = \J\uJ m \ 2 + 



v 2 F q 2 . 



(27) 
(28) 

(29) 
(30) 



An important feature of the kernel function K(iuj m 
uj m . Indeed, using S(iu m , q = 0) = W, So(iu m , 9 = 0) 



, q) , Eq. <|28(l , is that it is exactly zero for q = for arbitrary 
= |u> m |, and T(iu> m ,q = 0) = 1/(Wt- 1), we get 



K(iu> m , q = 0) 



1 

W 



1 



1 



Wt-1 
1 



1 

Wr 



1 



W-1/t \u m \ 



= 0. 



(31) 



We stress that this equality only holds when all the contributions to the self-energy are combined together. This 
property of the kernel function is characteristic for the gauge-invariant quantities in the presence of interaction. 
Indeed, the interaction at q = implies the shift of the chemical potential and hence can be gauged ouli^. Therefore 
the contribution of small q to the observables should be suppressed by vanishing of the corresponding kernel function. 
The same situation is well known for the interaction-induced correction to the conductivity^. 

It is worth discussing a peculiarity of the problem of magnetooscillations with respect to the gauge invariance. 
The gauge invariance of the oscillatory part of the observables is guaranteed by the fact that the thermodynamic 
potential is represented by closed loops. Since the characteristic spatial scale for such a loop is cyclotron radius R c , 
the interaction with momenta q <C R^ 1 should not contribute. So, if we would find that our result does not satisfy 
this requirement, it would mean that the diagrammatic treatment is not sufficient and should be complemented by 
the infrared cutoff at q ~ J?" 1 . On a more rigorous level this could be done by using the real-space path integral 
approach combining the treatment of magnetooscillations in the presence of long-range disorder— with the quantum 
kinetic-equation approach to the interaction effectsS^i and to magnetotransport 42 . In this context, it is instructive 



to recall the calculation of the dephasing length l v , where the infrared cutoff at q 



K 1 



or q 



R arises for the 



problems of weak localization^ and Aharonov-Bohm oscillations^ in a ring of radius R, respectively. 

Since we find, however, that the kernel K(iuj m , q) governing the perturbative self-energy l|27|l does satisfy K(iu) m , q = 
0) = 0, the above cutoff is irrelevant and the perturbative treatment is sufficient in the considered regime of strongly 
damped oscillations, T lo c . Indeed, the g-integral in l|27|) is cut off at q ~ lo/v f ~ T/vf which is much larger than 
R~ x under the above condition. 

We also emphasize that the kernel K(ioj m ,q) given by Eq. I|28|l vanishes in the clean limit 



KiitJ^q)^^ = 0. 



(32) 



This implies, in particular, that the correction to the effective mass found in Refs. I.32l8.'"l from the e- and p- 
dependence of the self-energy S(e,p) of a clean system, dm/m* ~ T/E F does not show up in the damping of 
magnetooscillations, in accordance with the statement made in Ref. .32. In general, S(e,p) is not an observable (and 
not gauge invariant) quantity, and thus the above correction 8m should be at least treated with caution. 
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III. DAMPING OF MAGNETOOSCILLATIONS: SHORT-RANGE INTERACTION 



In this section we evaluate (5E(i£ ra ,£o) i n the case of a weak point-like interaction given by 

V(iu m ,q) = U . 



(33) 



We are interested in the correction to the self-energy to the first order in uUq, where v = m* /2ir is the density of 
states per spin direction. We calculate the contributions E"- and starting from Eqs. i|B.15|) . i|B.8(l . (|B.9|) . (|B.20fl . 



(|B.21ll . and l|B.22|) derived in Appendix B. Using the notation S(iuj m ,q) introduced in Eq. (|29|l . the vertex factor T 
can be presented as 



1 



St-1 

Performing the momentum integration for < q < kp, we obtain 



(34) 



Sgo (»£»,&) = -iTU, 



2S^ (fe n ,&) 



qdq 
2tt 



1 1 

S Sq 



iTU x - 
2irv 2 ^ 

» u m >e n 

i T vUo <s~^ 
2E f t 2-" 

£?1 (^n,£o) = 



El 



W 2 ~W - \ E 



iTU 



iTuUo 
2E f t 



E 



2ir v F t 



J_ ,.,2 



,+l/r 



d5 



(2g~l/r) 
(5-1/r) 2 



(35) 



(36) 



These expressions are valid independently of the value of the parameter Tt, i.e. they describe both the diffusive and 
the ballistic regimes, as well as the crossover between them. The logarithmic term in l|36[) comes from Sg l5 while the 
l/u) m term originates from E^. 

The contributions of the Hikami-box diagrams are given by 



E b 



oo 



'ie n ) 



2E f t ^ 



S 2 



2E f t ^ 



- 1 



2ttv 2 t 2 ^ 



2E b Jl (ie n ) + T,\ 1 (ie n ) = 



u^+i S 2 (S-1/t) 



dS (2S-1/t) _iTvU a ^ 



2i?pT ' Ul m T 

<+>m>£n 



(37) 
(38) 



We see that the Hikami-box contribution exactly cancels the second (l/oj m ) term in 1)36(1 . Thus the total correction 
to the self-energy reads 



EpT 



- 1 



In 



(39) 



The upper limit of the summation over u> m is effectively given by m max ~ Ep/T. The term [uj m / y/E F + ui^ — 1] in 
Eq. (|39[1 . which originates from Eoo, yields after the summation over ui m > irT a contribution a T in addition to a 
large T-independent contribution, renormalizing r. As we will see below, a term oc T\n{E F /T) will arise from the 
contributions of Eoi and En, and hence Eoo yields only a subleading contribution to the damping. 
To calculate the sum of logarithms in (f 391) we write 



E * 



E 



l n ^ + l n V^ 



E F 



(40) 
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For T <C Ep the second sum can be replaced by the integral, yielding a T-independent contribution. We further use 
the identity 



M 



N 



El V 
In — = MlnN -hiT(M +1), 
m 



and the Stirling's formula 



hxT(M + 1) = M In M -M+~ ln(27rM) + ... 
for M 1, where T(x) is the gamma-function. As a result, we get for arbitrary positive e r . 



i T vU{) ( C\ Ep 



EpT 



T 



In — — — + In 



2ttT T 



1 _ r ( £r 



1 



/2tt \2ttT 2 



(41) 



where are constants of order unity. We see that the contribution to <5£(i7rT, £q) containing vertex corrections to 
the interaction line has an additional factor \n(Ep/T) as compared to the T-dependent part of (5Soo- ln Eq. i|41[l . 
we have absorbed the contribution of <5£oo m t° the upper cutoff of the log-term which is given by Ep up to a factor 
of order unity. Furthermore, the same can be done with the last term in Eq. i|41l which at e„ = irT also yields a 
linear-in-T contribution. Equation (|41|l thus translates for T lu c into the following expression for the damping 
exponent B(T) = — 2Tri5Y<(iiTT, (o)/w c : 



B{T) = -ci uU + - — In — • 

W C T Ul c hpT 1 



(42) 



The first term in Eq. (|42|l describes the T-independent FL-renormalization of r due to vertex corrections and should 
he included in the effective relaxation time r*. as was mentioned after Eq. 1151 in Sec.HTI The second term represents 
the T-dependent contribution to the damping factor that we are interested in and is analyzed in the next section. 



IV. INTERPRETATION: EFFECTIVE MASS VS QUANTUM SCATTERING TIME 



The above result (|42|l can be interpreted in terms of corrections to the effective mass (or u> c ) and the elastic scattering 
rate entering the standard formula (|17l) . These corrections come from the interplay of disorder and interaction. Writing 



■ exp 



47T 2 T 
UJc 

A { °\T) exp 
A ( °\T) exp 



27T 2 T 




IT 


1 


(uj c + 8uo c ) 


(u) c + 5u c ) (r 


+ St) 


2tt 2 T5u) c ~ 




TT f SuJ c 


St\ 


exp 






Ul c L0 C _ 




_U1 C T \ U! c 


'J 



2tt 2 T St 



m 



exp 



TT 
U1 C T 



8m 
m 



St 

T 



we conclude that 



B{T) 



2ir 2 T Sr 



n i 5m 
ui c t V m 



St 

T 



(43) 



(44) 



It is worth noting that the FL-renormalization does not affect the product uj c m = eB. 

Comparing Ij42(l and (|44|l [we recall that the first term in Eq. I|42|l is absorbed in r] , we see that the T In T dependence 
of the damping factor could in principle originate either from the InT correction to the effective mass, or from the 
T In T- type correction to r. This led the authors of Ref. l2a to the conclusion that the nonlinear T— dependence of the 
damping factor may be equivalently interpreted either as a T— dependent renormalization of the effective mass or as 
a T— dependent Dingle temperature. It is clear, however, that these two possibilities correspond to different physical 
processes. 
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A. Self-energy at real energies: analytical continuation 



To identify the physical origin of the leading contribution to the damping it is instructive to obtain B(T) using the 
expression for the self-energy analytically continued to real values of energies e n — > — ie . Performing the analytical 
continuation to real energies e and real frequencies w in Eq. (|27[1 . we get 



EpT 



doj 

47T 



tanh 



2T 



In 



(45) 



so that the real part of the self-energy is given by 



Re5E(e,£ ) 



AnEpT 



du> 



tanh 



tanh 



£ In 



e — uj 
2T 



2T 



tanh 



tanh 



e + u 



2T 



In 



2T 



hi 



2 - w 



yju 2 - E 



2-kE f t max[|e|,T]' 



(46) 



The leading contribution here comes from the term ]xl(Ef/uS) in the first integral over oj < Ep while other terms only 
rescale the ultraviolet cut-off Ep of the logarirthm by a constant of order unity, which is beyond the accuracy of our 
quasiclassical approximation. 
The imaginary part of ST, reads 



WE(e,£ ) 



AnEpT 



(IlU 



tanh 



2T 



- const - 



vU, 



ii 



T In 



r 2E f t 

where the T-dependent term has the following asymptotics: 

T ln[2cosh(e/2T)] = \ 



2 cosh 



tanh 

~) 

2TJ 



2T 



e/2, 
Tin 2. 



£> T, 
e<T. 



(47) 



(48) 



The contribution of the term u}/-*JE% — uj 2 to the integral in l|47fl is T- independent up to small corrections of order 
of vUq(T I Ep) 2 /t which are beyond the accuracy of the calculation. 

Having calculated Re S and Im E for real energies e, we can determine 8m and St. Indeed, the magnitude of the 
first harmonics of the magnctooscillations of the thermodynamic density of states is expressed through the real-e 
self-energy #E(e) as follows: 



MT) 



4n 2 T 



de 



dnp(e) 



Ai(e,T), 



2iri 



Ai(e,T) = exp<^ [e - ReS£(e,£o)] } exp <^ - 



2tt 

UJ C T 0J C 



Im<JE(e,&: 



(49) 
(50) 



where np(e) = [1 + exp(e/T)] _1 is the Fermi distribution function. 

In analogy with Eq. Ij43(l we represent the energy-dependent amplitude A\ (e, T) in terms of energy- and temperature- 
dependent corrections to the quantum scattering time and mass, 5t{e,T) and m(e,T): 



Ai(c,T) = exp 



2irie 



6m(e,T) 



exp 



7T 

LU C T 



Sm(e,T) 5r{e,T) 



Comparing 1)51(1 with l)50[l. we express 5r(e,T) and m(e,T) through ReE(e) and ImE(e) as follows: 

5m{e,T) Re,5E(£,T) 



m 
6r(e,T) 



2TlmSE(e,T) 



Sm(e, T) 



(51) 

(52) 
(53) 
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Using H46fl and (|47|l in combination with (|52|l and l|53|) . we obtain 



8m(e,T) uU Q ^ E F ( 



lirEpT max[|£|,T]' 



Sr(e,T) T 

= f£7 — In 

T Ef 



2 cosh f^)!--^ In E * v (55) 
V2T/J 2tt^t max[|e|,T] V y 



The integration in 149(1 sets in effect e ~ T in the above expressions. The T-dependent corrections to the effective 
mass and the quantum scattering time extracted experimentally with the help of Eq. I|43|) are thus given by 

(56) 



8m(T) 




1 Bp 
111 , 

T ' 


777, 


1-kEft 


Sr(T) 


T 

= Ep~ - 


vU Q 


T 


2nE F 



(57) 

In (|57f) we assume that ttvI/qT/Ef <C oj c t, expand exp[7r<5r(e, T)/uj c t 2 ], and then average the term ln[2 cosh(e/2T)] 
[which is a real-energy counterpart of the last term in Eq. (|41|l ] with ~dnp(s)/de. 

It is clear from these results that the leading term in B(T) [proportional to T\h(Ef/T), Eq. I|42|)] originates from the 
real part of the self-energy, Eq. JJBJ, i.e. from renormalization of the effective mass, which affects incommensurability 
of the oscillations at different values of energy e. The contribution to B(T) of the imaginary part of the self-energy 
[corresponding to the last term in Eq. 1(41(1 ]. which is governed in the ballistic regime by the renormalization of the 
scattering time, is smaller by a factor h\(Ep /T). In the expression for the damping, Eq. 142|l . this contribution is 
absorbed in the numerical constant Ci in the upper cutoff of the logarithm. 

The obtained result for the interaction- induced correction to the scattering time r, Eq. I|57[l . agrees, up to a 
factor i, with the correction to the transport time following from the calculation of conductivity correction in the 
ballistic regime in Ref. |£§ This is exactly what one would expect on physical grounds. Indeed, it is known that the 
conductivity correction^ can be understood as governed by an additional, predominantly back-scattering, contribution 
to the scattering cross-section related to the dressing of an impurity by Friedel oscillations. Since this contribution is 
concentrated near the scattering angle <f> = tt, the correction to the momentum relaxation rate is larger by the factor 
1 — cos0 ~ 2 than the correction to the total scattering rate. In Sec. HVBl we will corroborate the results of this 
subsection by an explicit calculation of the contribution to the impurity scattering rate due to Friedel oscillations. 

Up to now we calculated the exchange contribution to the self-energy. For the point-like interaction, the Hartree 
term has opposite sign and is twice larger in magnitude than the exchange term due to the spin summation. This 
simply reverses the sign of the corrections to the damping factor. 



B. Calculation of St from the scattering off Friedel oscillations. 

In this subsection we calculate the correction to the total elastic scattering time 8t in a different, physically more 
transparent way, considering the scattering off impurities dressed by Friedel oscillations^^. We will demonstrate how 
the result (|57|l is reproduced in this way. In particular, this will confirm once more that there is no TlnT term in 8t 
and therefore the leading Tin T contribution to the damping factor comes from 8m. 

We start with the expression relating the total elastic scattering rate and the scattering cross-section S{(j)) of a 
single impurity, 

vf I P-Stf), (58) 



r q (e,T) * r J 2n 

where nj mp is the concentration of impurities. The expression for the transport scattering time Tt r determining the 
conductivity differs from l|58l) by a factor 1 — cos <f> in the integrand: 



1 



rtr{e,T) 



/f;<SW)(l-oos0. (59) 



We note that the two times, r q and r tr , though equal for the point-like impurities in the non-interacting case [for 
which S(<fi) — (jT.impTjpT) -1 = const(0)], differ from each other when the scattering off Friedel oscillations is taken 
into account. 

The impurity scattering cross-section for dressed impurities reads 



s(<t>) = 2 — 

VF 



V ( 2k sin | 



~ S + -^!-JL 2 Re <5U ( 2k sin - I . (60) 
v F V 2 
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FIG. 4: Schematic illustration of the scattering off Friedel oscillations. The black dot in the middle represents a short-range 
impurity which creates the oscillatory correction to the electron density around it. The circle represents the equipotential line 
of the effective impurity potential. The correction to the impurity cross-section at the angle cj> arises due to the interference of 
two electronic waves*!*— , one of which (dashed line) scattered by the impurity and another (solid line) by the Friedel oscillations 
at a point parametrized by the distance r from the impurity and the angle %j>. 



Here So = 2hvVq/vf is the bare impurity scattering cross-section and Vq = J d 2 rVo(r), where Vb(r) is the bare 
point-like impurity potential. The cross-section S((f>) depends on energy e of an electron through k = kp + e/vp in 
the Fourier transform V(q) = Vq + SV(q) of the effective impurity potential, renormalized by the Friedel oscillations 
of the electron density. For r ^> kp 1 the oscillatory correction to the electron density reads 



Sp(r) = -vV — 



(2TrrT/vp) sm2kpr 



smh(2TrrT / v f) irr 2 



(61) 



The correction to the impurity scattering potential SV(r) due to scattering off the Friedel oscillations is proportional, 
for the short-range interaction, to the electron density at point r. Similarly to the consideration of Sec. IIIII and 
Sec. IIV Al we will concentrate on the exchange part of this correction, 



5V{r) 



^U Sp(r) 



To calculate the correction to the impurity cross-section, we need the Fourier transform of 5V(r), 

SV ^2fcsin^ = ^ J d 2 rSp(r) exp ^2i krsin — cos'ipj , 



(62) 



(63) 



where ip is the polar angle of r, see Fig. 0] Substituting (f^T~f) and ljtj5|) into (|fJU|l . we find the interaction-induced 
correction to the scattering cross-section, 



5S(<t>) = uU S 



2^ 



dr- 



(2ttT/v f ) 



sinh(2nrT/vF) 
Performing the integration over r we get 



sin 2r 



kp + k cos ip sin ■ 



sin 2r 



kp — k cos ip sin ■ 



<JS(<£) = S V —£- J d^j [tanh (3+ + tanh/3_] , 



where 



p± = — kp ± k cos tp sin — = — — 



2T 



T 



1± 1 



2E F 



cos %b sin — 
2 



(64) 



(65) 



(66) 



For j3± 3> 1, i.e. for most values of the scattering angle 4> except for those corresponding to the backscattering 
((j> fts 7r), we see that the scattering cross-section does not depend on ef> up to exponentially small corrections of order 
0(exp[-E F /T\): 



5(0) = So [1 + nisUo] , \<j> - ?r| > (max[T,e]/E F ) 1 / 2 . 



(67) 
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FIG. 5: Differential impurity cross-section S(4>) (in units of the bare cross-section So) renormalized by the scattering off Friedel 
oscillations (with only the exchange contribution taken into account) calculated for uUo = 0.2 and (a) fixed T = 0.01-Ef and 
e/E F = 0, 0.01, 0.02, 0.05, 0.1, 0.2, from top to bottom; (b) fixed e = and T/E F = 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, from 
top to bottom. 



On the other hand, in the vicinity of = 7r the cross-section has a "hump" , see Fig. \5\ of the width and the height 
scaling as <50, 8S(4>)/S ~ (max[T, e]/E F )^ 2 . The explicit expression for S(4>) at T — can be found in Ref. U In 
Fig. we plotted S(cj>) for fixed T and several values of e (Fig. Et) an d for fixed e = and several values of T (Fig.Eb). 

To calculate the (e, T)-dependent correction to l/r q (which is determined by the "hump"), we expand sin(</>/2) 
around <j) = n. Furthermore, we expand cost/' around ip = tt and ip = in the expressions for f3 + and /3_, respectively. 
This corresponds to the two possible interfering paths propagating along the horizontal line in Fig.^lwith the scattering 
off the Friedel oscillation occurring either to the left or to the right from the impurity. Denoting Sip = x, 8<p/2 = y, 



y 2 , and u> = z 2 E F , we get from Eq. I|65[l 



2tt 



SS^) ~ So 



vUo 



2tt 
S vUo 



dx / dy tanh 



T 



2E F 



2 i 2 

x +y 



zdz tanh 



2E F 



So 



vUo 
2E F 



duj tanh 



2T 



Using l|58|l and (|68|l . we find the total quantum scattering rate, 



Ml 



«T,(e,T) 



r 



const In (2 cosh — ^ 

E F V 2TJ 



2E F 



(68) 

(69) 
(70) 



This result for the interaction- induced correction to r q agrees (to the leading order in 1/Tr corresponding to the 
ballistic regime) with that obtained from the imaginary part of the self-energy Eqs. (|55ll and (|57|l . More accurately, 
Eq. (|70|) differs from Eq. I|55|l by the last (T-independent) term vUo£/2E f t, which drops out after the thermal 
averaging with —dn F (e)/de and thus does not contribute to Sr q (T), Eq. I)57p. This term is in a sense anomalous, 
since it arises from the ultraviolet limit of the cj-integration in Eq. I|68() . In fact, one could question the validity of 
this contribution, since we used the asymptotic, large-r form of the Friedel oscillations in Eq. I|61|) . One can check, 
however, that the same result [up to an irrelevant additive constant independent of T and e: in Eq. (|67J) itvUo is 
replaced by 2vUo\ is obtained from a calculation using the exact form of the Friedel oscillations^. The appearance of 
this linear-in-e term is related to the violation of the particle-hole symmetry in the parabolic spectrum; this term did 
not appear in the diagrammatic calculation of Sec. II VI where the spectrum was linearized. What however enters the 
experimental damping of magnetooscillations is l/r q integrated over the energy with an even function —dn F (e)/de. 
Therefore, we are in fact interested in the collision rate symmetrized with respect to e — * — e. Performing this 
symmctrization in Eq. (|64|l . we get a result determined solely by the infrared scale, r ~ v F /T, yielding Eq. I|70() 
without the last, linear-in-e term. 

Finally, using l|59|) . we see that the correction to the transport rate is larger than <5(1/t 9 ) by a factor of 2. The 
corresponding correction to the conductivity reproduces the result of Ref. .9, in the ballistic limit. 
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FIG. 6: Diagrammatic equation for the effective interaction line (bold wavy line) in the random-phase approximation, Eqs. I|71|l 
and 1731 1. Dotted wavy line represents the bare interaction, Vo(q) + F ( p jlv in the singlet channel or F§ /2v in the triplet channel. 
The bubble II is the polarization operator, Eq. I72H . 



V. DAMPING OF MAGNETOOSCILLATIONS: COULOMB INTERACTION 



We turn now to the Coulomb interaction. In the case of Coulomb interaction, one should take into account the 
dynamical screening of the interaction within the random phase approximation (RPA), see Fig. El In what follows 
we use for simplicity the so-called i^-approximation^, which retains only the zeroth harmonics of the Fermi-liquid 
constants Ffi and Ffi in the charge and spin channels, respectively. Then the effective interaction propagator in the 
charge channel (combining the exchange term and the singlet contribution of the Hartree-type interaction) readsSiIfl 

V"(iu m ,q)= [(V (q)+FP/2u)- 1 +Il(iuj m ,q)]- 1 , (71) 

where II (iu) m , q) is the polarization operator 

II (iu m , q) = 2v [1 - |w m | tT {ioj m , q)] , (72) 

r (iuj m ,q) is the impurity ladder, and v = m/27r is the density of states per spin direction. The triplet contribution 
to the effective interaction arises from the ladder of Hartree-type interaction blocks and reads 

V a {iu m , q) = [2v/F° + n (tw mj qT 1 . (73) 



A. Singlet channel 

The main difference as compared to the case of the short-range weak interaction considered in Sec. Ill is the 
nontrivial form of the dynamically screened Coulomb interaction (one should take into account the renormalization 
of the interaction by polarization operator) . Using 

Vo(q) = 

q 

and neglecting F°, we get 

V (q) 2we 2 /q 



V p {lLO mi q) = - 

J- "1" KO^MUmji/j ± t yz,-/it- / q)z,u[± — \w m \-i i yiu/ m ,q)\ 

(74) 



2v (q + k)(S - 1/t) - k(W - 1/r) ' 

Here we use the standard notation k = Airve 2 for the inverse Thomas-Fermi screening radius and use the short-hand 
notation S for S(iu) m , q) = \JW 2 + v 2 F q 2 with W = \uj m \ + 1/t. 

For q <C K, neglecting q in the sum q + k in the denominator of l|74|l. one finds that the exchange interaction 

y p (^) = h^w (75) 

has a singularity oc l/q 2 in the limit q —> 0: 

V p (tuj mi q ) = - — , q^O, (76) 

2v q Vp 

so that each separate term £■"-, in the self-energy would diverge. This divergence is analogous to the one 

encountered in course of calculation of the tunneling density of stateai In that case, one has to keep the q term 



15 



in denominator of 1|74|) . which cuts off the logarithmic divergence. For the present problem, this is, however, not 
needed. Indeed, as was emphasized in the end of Sec. ^ the kernel K(iu> m ,q) combining together contributions 
of all relevant self-energy diagrams is proportional to q 2 in the limit q — > and hence cancels the singularity in 
V(ito m , q). In view of this, it is convenient to represent the kernel function H28|) in a form which shows explicitly that 
K(iu m , q = 0) = [we recall that S a (q = 0) = \cj m \ and S(q = 0) = W}: 



K{iu mi q) = [ -i + — !- 



Then the product V p (iui m , q)K(itu m , q) takes the form 



S-W 
S-1/t 



1 



1 



2vV p (iuj m , q)K(w m , q) 



So 



St(S-1/t) \u m \ 



1 



So(S-W) S St(S-1/t) 



(77) 



(78) 



Performing the integration over the momentum q in i|27|) , we get the correction to self-energy in the singlet (charge) 
channel 



i T 
2E f t 



A r 



2uj m T) In — — — 1 + In 



2uj m T 



where we introduced 



KVf- 



(79) 



(80) 



Comparing (|79|l and l|39[l . we see that the last term in (|79() corresponds to a static short-range interaction with 
uU = l. 

Setting e n — -kT and separating the contributions to the sum l|79|l governed by the high-energy (e ~ A) and 
low-energy (e ~ T) regions, we can present the result in the following form: 



where f(x) is a parameterless function, 



i T 
2E f t 



constA 



T 



1 



1 



8itTt 



hi- 



/(*) = E 



m— 1 



1 — (1 + mi) In 



1 + mx 



2mx 



Cf 

x 



T 
1 

2x " 



/(47rTr) 



(81) 



In- 



36 x 2 ' 



x < 1, 



x > 1, 



(82) 



with Cf = —3/4 — ip(l)/2 = —0.461392.. (here ip(x) is the digamma function). Thus the dynamical screening of 
Coulomb interaction leads to different asymptotics of the self-energy in the diffusive and ballistic regimes, in contrast 
to the case of weak short-range interaction^. 

The T-dependence of the leading correction to the magnetooscillations damping factor due to the interaction in the 
singlet channel has therefore the form 



B p (T) 



7T T 
U) C T Ef 



1 



8ttTt 



ln- + /(47rTr) 



(83) 



, -In— - - ( H \ — ] ln(47rAr) —, AttTt < 1. 

— — x <^ (84) 
uj c t E F / i \ A 

1 In—, 4ttTt>1. 

8itTt J T 

This result is illustrated in Fig. Note that the ballistic asymptotics describes the exact result with a remarkable 
accuracy down to very low temperature, Tr ~ 0.01 — 0.05, see Fig. 01. Retaining in Eq. i|84J) only the leading terms 
and suppressing the T-independent contributions which can be absorbed in the FL-renormalized r, we get 

f ^ ln ^-^ ln (47rAr), AttTt « 1, 

B»(T) = — — x{ ~ (85) 

lo c t Ef I A 

[in—, 4ttTt>1. 



16 



a 
P 

CL 

CP 




0.5 



-i 

e 

w 

3 
P 

CL 

PQ 



-0.5 



<* / 
/ 



Tx 



0.1 



FIG. 7: Temperature dependence of the singlet channel correction to the damping factor B P (T), Eq. 1831 . for 4ttAt = 100 
(solid line) with the low-T (dot-dashed) and high-T (dashed) asymptotics, Eq. 18411 . (a) Wide temperature range: on this scale 
B P (T) is essentially indistinguishable from its high-T asymptotics; (b) low-T part: the crossover between the two asymptotics 
occurs at Tr ~ 0.05. 

B. Triplet channel 

Calculation of the corresponding triplet contribution B a (T) is presented in Appendix C and leads to qualitatively 
similar asymptotics. The leading term in the total correction to the damping factor in the ballistic regime, realized 
in experiments on low-disorder samples at realistic temperatures, takes the simple form 



B(T) = B P (T) + B a (T) ~ 1 



- + Fo 



7T T A 

oj c t Ef T 



(86) 



As discussed in Sec. II VI this result arises due to the correction to the effective mass (the consideration of Sec. IIVI fully 
applies to the case of the Coulomb interaction as well). 

It is worth noting that due to the factor (1 + Ffi) in the denominator of B a (T), the damping of magnetooscillations 
tends to diverge upon approaching the Stoner ferromagnetic instability. Since the damping is determined by the 
effective mass m*(T) (see Sec. IIV|) . we conclude that m*(T) diverges as Ffi — ► —1 due to the interplay of disorder 
and interaction. This should be contrasted with the clean case, where the effective mass is solely determined by 
Fermi- liquid constant in the singlet channel, Ff, and hence is insensitive to the ferromagnetic instability. 



C. Discussion and comparison to earlier work 



Let us discuss the obtained results for the damping factor B(T). In both the diffusive and ballistic regimes we find 
the temperature dependence of the form B(T) oc TlnT. In the diffusive (low-T) regime the relative correction to the 
damping factor, u> c tB ~ g _1 Tr In T is less singular than the known corrections^ to the conductivity, 8a j a ~ <? _1 InT, 
and the tunneling density of states, 5p/p ~ <7 _1 ln T (here g ~ EpT is the dimensionless conductivity). On the 
technical level, the qualitative difference with the tunneling DOS can be traced back to the contribution of the 
Hikami box diagrams, Sy. 

As has been already mentioned in Sec. [I] the problem of the effect of the interaction on magnetooscilations was 
recently addressed in Ref.|26|. The result of this work for the damping factor is qualitatively similar to ours, B oc T In T. 
However, the crossover function and, in particular, the prefactors in both ballistic and diffusive limits differ from ours 
(1 instead of 3/2 in the first line of Eq. I|85|). and 3/2 instead of 1 in the second line). This difference is not surprising, 
since the authors of Ref. |26| took into account only one diagram out of six in Fig. [2] Thus, even the qualitative 
agreement may be considered as an accidental coincidence. In fact, there is a conceptual difference between our result 
and that of Ref. [2(| To illustrate this, consider a toy interaction of the type V(u>,q) = F(uj)S(q). Our result then 
would be zero, since the kernel function K(iuj m ,q), Eq. satisfies the gauge-invariance constraint [see Eq. (fUTf) 
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and discussion below it], 

J K'(iw ro ,q)|q=o = 0. (87) 
In contrast, the formula of Ref. would give a finite result, since their kernel 

K {a) --- 1 

-"■01 — c 



2 



SSt- 1 

does not satisfy the requirement l|87|) . 

K$(iu m ,q= 0) = — 

\u m \{\Um,\T + 1) 

For the Coulomb interaction this results in a logarithmic divergency at small q that is cut off by the plasmon pole at 
q ~ q min with 

_ |tJ m |(|a; m | + 1/t) 
KVjr 

similarly to the calculation of the tunneling density of states^. As we explained in the end of Sec. [HI the contribution 
of small momenta, q < R ~ , should be suppressed for the present problem. Therefore, the above small-g divergence 
should be cut off by the magnetic field, which would partly transform InT of Ref. [2|3 into ln_B. However, this problem 
is in fact spurious: the result of our work does not suffer from any infrared divergencies, since our kernel K(iuj m ,q) 
does satisfy Eq. (|8*7jl . 

Finally, let us briefly comment on the ultra-diffusive regime, T <C uj c <C 1/t. In this regime the summation over 
Matsubara energies e n is not restricted to n = and N c ~ lu c /T 3> 1 Matsubara harmonics are important. Therefore, 
the damping due to the inelastic scattering [suppressed only for n = 0, see Ref. and Eq. (|B. 13|l ] becomes finite. 
The corresponding contribution to the damping can be roughly estimated using Eq. (|B.13(I taken at relevant n ~ N c : 



T 2 

SB incl (T) ~ - n 2 \u(E f t) 

LU c Ep 



~-^Lln(E F r), (89) 



yielding \5B lncl (T)\ <C 1, since in the ultra-diffusive regime w c < 1/r. Thus the inelastic contribution to the damping 
factor is always small. Note that the contribution to the damping due to the renormalization of the effective mass in 
the ultra-diffusive regime is also small: at n = N c we have B(T) = —2ite n /u) c 6m/m ~ Sm/m <C 1. 



VI. CONCLUSIONS 



In conclusion, we have studied the T-dependent interaction corrections to the damping of magnetooscillations in 
a two-dimensional electron gas. The damping factor has been calculated for Coulomb and short-range interaction in 
the whole range of temperatures, from the ballistic to the diffusive regime. While the relevant diagrams are similar 
to those for the local density of states, the results are essentially different, see Eqs. 142(1. JH3J|. 

We have identified leading contributions to the damping induced by interplay of interaction and disorder, which can 
be associated with corrections to the effective mass and the quantum scattering time. It has been shown that in the 
ballistic regime, which is typically realized in low-disorder samples at realistic temperatures, the dominant effect is that 
of the renormalization of the effective electron mass due to the interplay of the interaction and impurity scattering. 
Specifically, the correction to the effective mass is of the form 5m/m ~ 1/(Eft) ln(Ep/T), Ea. (|56|l . The correction 
to the impurity scattering time is of the form Sr q /T ~ T/Ep, Eq. Ij57[l. and yields a subleading contribution to the 
damping. We have confirmed the result for the correction to the quantum scattering time by performing a calculation 
based on the picture of scattering by impurities dressed by Friedel oscillations. The results of the paper are relevant 
to the analysis of experiments on magnetooscillations (in particular, for extracting the value of the effective mass) and 
are expected to be useful for understanding the physics of a high-mobility 2DEG near the apparent metal-insulator 
transition. 
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APPENDIX A: LUTTINGER FORMALISM FOR THE THERMODYNAMIC POTENTIAL 

In this Appendix we derive the oscillatory part of the thermodynamic potential 0, following Ref. l3~il 
We calculate the sum in Eq. 0, using Poisson's summation formula 

f(^c) = 5/(0) + / ^ f(0 - E o-T / # f'(0 (e^M - e- 2 -^/-) . (A.l) 

m=0 ^° Wc fe=l 71-1 ^° 

From Eq. (J7J) we have /'(£) = — G(ie„,£), where G(ie„,£) is the Green's function. Thus, to extract the oscillatory 
contribution to f2, we should calculate the following integral 



Jo £ - M - + £ co (ie„,£) - ?sgne„/2r [ <9£ 



(A.2) 



Here we introduce the self-energy Y, ee (ie n ,£,), which is a function defined in the plane of a complex variable £, such 
that 

Ecc(«£n, mWc) = £cc(«£n, £ = ™w c ). (A. 3) 

The main contribution comes from the pole £ = £bj where £o obeys the self-consistent equation 

i 

£o = H + ie n + — sgne„ - E co (ie„, £ )- (A. 4) 

We expand the self-energy (from now on we will skip the subscript "ee" in E C c) in the vicinity of the pole, 

<9£(i£„,£) 



9£ 



(A.5) 



Then the denominator in Eq. 1A.2|I becomes proportional to (1 + 9E/9£) so that these factors drop out and the 
integral takes a simple form 

J a ? — so 

We first single out the FL renormalization factors in X, i.e. represent the self-energy in the following form (assuming 
a constant electron concentration) 

E(ie„, £) ~ Sfx + /3 (£ - p) - ia e n + ST,(ie n , £), (A. 7) 

where 

fi = fj, — SfJ, = Tm e /m 

is the chemical potential for noninteracting electrons (n e is the electron concentration): £ — p ~ Vp(k — kp). The 
correction to the self-energy ST,(ie n , £) contains contributions that are smaller that /3o(£ — /•*) and ao£n by either T/Ep 
or 1/ Eft. These additional contributions are related to the inelastic processes and to the modification of the pure 
FL result due to disorder. 

In order to solve Eq. (|A.4|I . we treat the subleading terms constituting as small corrections. Solving Eq. I|A.4|) 
by iterations, we first find its solution neglecting these small corrections and expanding Eqq around the mass-shell: 



which yields 



?o 0) - A + i£ n + 7T sgne„ + ia e n - /3 (£o 0) - p), (A. 
At 



AO) _ Th. + j F 1 + a ° + lSgIi£n (A q) 

Co _ m +Ze "l + /3o + 2r(l + /3 )- (A ' 9) 



Next we use this value of £o i n and solve the self-consistent equation again, now keeping the terms previously 
neglected. Then we arrive at 

£ ~ - 5E (ie n ,^) /(l + [3 ). (A.10) 
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Substituting this value of the pole in Eq. (|A.6|) . we obtain 

^ — ^osc ^Non- osc i 

where 



(A.ll) 



and 



n— — oo 



Tvuj c In ( -/i - ie„ + S(ie„, 0) 



zsgn£ r . 
2r 



n— — oo 



2vT / dein(C-M-fen + S(*en.O- 



asgne n 
2r 



OO OO 



EV- T2vu c (2nkf.im e |e„|(l + a ) + 1/2t + iSE(ie n , f ) sgne„ 
V — r — exp<^ U sgn £ „ 
^-^ K W f \ 777, 
fc=ln=-oo 



fe=l 



k~ ° 0S eB ^ eXP \ uu c 1+00 

w c \ 2 27r 2 n e 
2/y — Ai cos — — — . 
\2tvJ eB 



1 + A> 
1 7<5S(ze n ,^ ) 



2r(l + a ) 1 + a 



Here A\ is the amplitude of the principal harmonics of the oscillations, 



4tt 2 T 



± >r-^ / 2ir 1 + a 



1 



^S(i£n,Co) 

2T(l + a ) l + a 



(A.12) 



(A.13) 
(A.14) 

(A.15) 



APPENDIX B: CALCULATION OF SELF-ENERGIES 

In this Appendix we calculate the relevant self-energy contributions and derive Eq. (|28Jl . The zero-B Green's 
function is given by 



G(ie n - iu m ,p- q) = 



{uj m - e n )(l + £*o) + sgn(w m - £„)/2t - - p) - iv F qcos6](l + /3 ) 



u) m -e n + sgn(w m - e„)/2r - i(£ - p,)(l + (3 ) + iv F qcos0 



(B.l) 



We denote the FL-renormalized energies and momenta as u> m = (1 + ao)w m , £ n = (1 + ao)e n , and q = (1 + /?o)<Z- As 
the first approximation, we have set in Eq. (|B.1|) p/m = ?jp in the linear-in-g term ivFqcos9 and neglected q 2 /2m. 

Since the effective interaction (I71|) and (|73|) does not depend on the polar angle of the transferred momentum 
q, we average the FL-dressed Green's function G(ie n — ioj m ,p — q) over the angle between p and q. The result of 
angle-averaging is 



d6 

(G(ie„ — iu m ,p- q)) = I —G(ie n - iu m ,f— q) 

Zn 



i sgn(w m - e n ) 



Vpro - e n + sgn(w m - £„)/2r - i(£ - p,)(l + O )} 2 + v 2 F q 2 



(B.2) 
(B.3) 



We substitute £ = £ ^ £o° } = M + + A)) + «sgn£„/2r(l + /3 ) for £ in since we are interested in 

<5£(i£ n ,£o) (the only place where £o appears in <5£ is the Green's function under the interaction line). Then the 
denominator in l|B.3(l for q = reads 



£o m -e n + sgn(w m - £„)/2r - i(£ - p,)(l + (3 ) 



u m -i n + sgn(w m - £„)/2r - i(l + O ) 



sgn £„ 



1 + 00 2r(l + /3 ) 



ti m + ~0(uj m - e n )6(e n ), 

T 



(B.4) 
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with 9(x) the theta-function. For definiteness, below we consider e n > 0. 

For S 01 and Sn at e n > we consider the Green's function at oj m — s n > in order to have different signs of 
Matsubara energies in the Green's functions connected by the interaction vertex. This condition allows us to dress 
the interaction vertices by impurity ladders. 

We see that e n drops out in the averaged Green's function taken at £ = £oi as in two-particle quantities: 

(G(ie n - iuj m ,p- q}} | £=£ = == , uj m > e n . (B.5) 

V(l w ">l + 1 / T ) +v 2 F q 2 

When both e n and e n — oj m have the same sign (such a contribution appears in the calculation of (5Soo), we find 

i 

yj\ib m \ 2 + v 2 F q 

Now we re-define the Fermi velocity to absorb the FL-factors according to 



(G(ie„ - iuj m ,p- q)) | £=£o = y=====, u> m < e n . (B.6) 



„ 1 + A) k F 

1 + ao m* 

Then we can return from uj m and q to uj m and q, expressing the angle-averaged Green's function in terms of FL- 
renormalized parameters Z , r* [introduced in Eq. i|15|) ] and v* F : 

iZ 

(G(ie n - iw m ,p- q)} L £ = = , u m > e„, 

4 40 \AHr1 + 1/r*) 2 + (vj^) 2 

(G(fe„ - iw m ,p- q)) | £=£ = ==, w m <£„. 

Furthermore, the Z-factor will be cancelled in the final result, when (G(ie n — iu> m ,p — q)) is used to calculate the 
correction to the observables, see e.g. Ref. 39 and discussion in Sec. IIB. 
Using (|B.5|) and i|B.6(l . we obtain 

/d 2 q 
—2 V(iuj m , q)T(iw m , q) (G(ie n ,p - q)) L £o 
, ( 27r ) 

cPg F(iw m , q) T(iu m , q) 

/ / 



^ / (2tt) 2 S-(ic; mj q) ' 1 



s?i(^„,Co) = -t [ ^^K.q)r 2 K,q) <G(^„,p-q)>| £ _ £o 



W m >£ n ' 



(2^ 

d 2 g U(iw m ,q)r 2 (iw m ,q) 



where for brevity we introduce new variables S(iuj m , q) = ^/ (|w m | + 1/t) 2 + v F q 2 = \JW 2 + v F q 2 and W = \u m \ + 
1/r. 

Now we consider the contribution to the self-energy without vertex corrections, £g (i£„,£). We recall that in 
£[)o(* £ n,£) ^ ne summation over transferred frequencies is not restricted to u m > e n . Presenting Sg (ie„,^) as 

SSo(*en,eo) = SSo +_ (^^o) + S^ ++ (z£„,Co), (B.10) 

/d 2 q 
——2 V(iw m , q) (G(ie n , p - q)) L £ 

T ^ / (^f%^' (B - n) 



-2 



we further split the contribution £qq ++ (ie n , £q) corresponding to no change of Matsubara frequencies at the interaction 



vertices into two parts as follows 



Zoo ++ (ie n ,£, ) = -T ^2 7^2 V ( iuj m, q) (G(ie n , p - q))| £ £o 



^T [**\ Y + Y Ittd (B.12) 

J ( 2 -) L,fe„ « m <-J V\"m\ 2 + q 2 v 2 F 
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The terms with \ui m \ < e n are responsible for the FL-renormalization and for the inelastic (determined by real 
processes) contribution to the self-energy, yielding the following FL-type term2&: 



E^ten^o) 



-ia 



i-y(ie n ,T)- 



ir 2 T 2 



(B.13) 



where the function ~f(ie ni T) depends logarithmically on max[e„, T, 1/r]: in particular, 7(ie„,T) oc \n[Ep/(— ie n )] for 
e n ^> T,1/t and 7(ie„,T) cx ln[ii?i?r] for 1/t 3> T,e n . The first term in IB. 13(1 determines the FL Z-factor and 
has been separated from 5E which governs the correction to the damping factor, see Eq. (|A.7|) . As for the second 
term, its imaginary part describes the inelastic electron-electron scattering, while its real part contributes to the 
renormalization of the effective mass^Si^. However, when taken at eo = ?rT, as appropriate for the damping of the 
magnetooscillations at T ^ u> c , the second term in l|B.13|) vanishes, in agreement with Ref.|2(J Note that for the case 
of weak short-range interaction, V(iu> m , q) = const(w m ), the inelastic contribution is zero to the first order in V. 

Thus in order to calculate <5E, we shall retain in (|B.12I) only the term corresponding to the summation over 
w m < —£ n . In this term we change the sign of io m and (suppressing the irrelevant inelastic term) and obtain 



d 2 q V(iu) m ,q) 



q 2 v\ 



thus arriving at 



1 



S(icj m ,<l) S (iu> m ,q) 



(B.14) 



(B.15) 



Here we use the fact that V(iuj mi q) — V(— iw m ,q) as the dynamically screened interactions depends on tv m only 
through \u! m \ and introduced So(iuj m , q) = y 'ui^ + q 2 v 2 F . 

Let us turn now to the Hikami-box diagrams, shown in Fig. [2t> We remind the reader that these diagrams are 
generated by covering each contribution to the self-energy Ey from Fig. [2K with e n (e n — tu m ) < by a single impurity 
line. Therefore for white-noise disorder (addressed in this paper) the Hikami-box contribution to the self-energy is 
independent of £ p and can be expressed through the corresponding E?- as 



<% p v[G(ie n ,Q] 2 ZUie n ,Z P ) 



i 



(B.16) 
(B.17) 



i d 

since the pole of G is given by £0, see Eq. (|A.10(I . Differentiating (|B.3|) we get 



This yields 



— G(ie n - ioj n ,p- q) 



Eoo(* £ "> Co) 
Eoi(fenj £o) 



T E 

<*> m >e t 

T E 



|a; m | + l/r _W 
{(\cu m \+l/r) 2 + v 2 F q 2 } 3 / 2 



d 2 q V(iu m ,q)W 



(2tt) 2 rS 3 {iuj m ,q) 
d 2 q V(iuj m ,q}T(iuj m ,q)W 



(2tt) 2 rS 3 {iuj m ,q) 
V 6 /■ tf N ■ ™ /" A y(^ m ,q )r 2 (^ m ,q) W 



^(iw^q) 



(B.18) 
(B.19) 

(B.20) 
(B.21) 
(B.22) 



Combining all the contributions together, we arrive at Eq. (|28|) . 



22 



APPENDIX C: COULOMB INTERACTION: TRIPLET CHANNEL 



In this Appendix we calculate the contribution of the triplet channel to the damping of magnetooscillations. The 
effective interaction in the triplet channel can be found by replacing Vo(g) — ► Fq /2u in the expression for the singlet 
channel: 



2vV a (ioj m ,q) 



1 + F$Il(iu m , q) 1 + F£[l - \u m \/{S - 1/t)] 

fs(s-i/t) 



(1 + F$)(S - 1/t) - \iu m \F° 1 ~i~ Fq S-w ' 



where we introduced 



w = \ui m \ 



1 + Fg r 



This yields 



2v V a (iu m , q)K(iw m , q) 



1 + Fg 



So -(S- 1/t) S-W 



Sq(S — w) S — w St(S — 1/t) 



(C.l) 



(C.2) 



(C.3) 



Performing the integration over q in l|27|l and taking into account the three triplet terms corresponding to different 
projection of the total spin on the z-axis S z = 0,±1, we obtain 



5S ff (i7rT,Co) 
h(z,y) 



i T 3F CT 
2E f t 1 + F ff 



E 



l =2irT 



- 1 



1--) in[i + y] + _X- zht{x,y) 
VJ 1 + 1/ 



In- 



+ /i(w m T,F ( f)|, (C.4) 
(C.5) 



h\{z,y) = In 



(l + 2z)(l + y) 
2z 



zy+(l + y) 



. zy-1 . zy+(l + y) 

arcsm , — arcsm ■ 



y/2z + 1 



X < 



In 1 



(l + y)V2z+T 
2 + 2/ \ 



zy - 1 + v /(zy) 2 -2(l + y)z 



- ln(l + y) 



z < 



z > 



2(1 + 1/) 
V 

2(1 + 1/) 

y 2 



(C.6) 



We see that the first two terms in (IC.4II correspond to the point-like interaction with i/Uq — * 3Fq /(I + Fq), see 
Eq. i|39|) . The term h(u> m , Fq) corresponds to the crossover function f(x) in the singlet channel, see Eqs. i|81|) and 
(|82|l . The result for the singlet channel is reproduced in the limit Fq — > oo (cf. Ref. 0j). 
The summation over Matsubara frequencies leads to 



<f£ <T (i7rT,£ ) 



T 3^ f constA 



2^ti + f ct \ r 



A(F CT ) 



8ttTt 



ln^-/ CT (47rTT,F CT ; 



where 



A(l/) 



1 + 1/ 



[y(6 + y)-2(3 + 2y)ln(l + y)] 



(C.7) 



(C.8) 



and 



OO s 

m=l ^ 



1 - 



21n(l + y) 



MX 
2mx J 



(C.9) 
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As a result we obtain the following T dependence of the triplet contribution to the damping exponent 



B"(T) 



7T 3F CT T 



8ttTt 



In- 



T 



U(4ttTt,FZ) 



(CIO) 



For not too strong interaction, (Fq ) 2 /(l + Fq) < 1, (i.e. for \Fq \ < 0.6, which is typically met in experiments, see 
e.g. R.ef. Efih the crossover function f a (x,y) only yields the subleasing T-dependence^ of B a (T), so that the leading 
contribution to the damping is given by Eq. ljl2"|) for the short-range interaction with vU$ replaced by 3Ffi /(l + Fq). 
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